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Abstract 

We perform a detailed study of the effective Lagrangian for the Goldstone mode of a superfluid 
Fermi gas at zero temperature in the whole BCS-BEC crossover. By using a derivative expansion 
of the response functions, we derive the most general form of this Lagrangian at the next to 
leading order in the momentum expansion in terms of four coefficient functions. This involves 
the elimination of all the higher order time derivatives by careful use of the leading order field 
equations. In the infinite scattering length limit where conformal invariance is realized, we show 
that the effective Lagrangian must contain an unnoticed invariant combination of higher spatial 
gradients of the Goldstone mode, while explicit couplings to spatial gradients of the trapping 
potential are absent. Across the whole crossover, we determine all the coefficient functions at the 
one-loop level, taking into account the dependence of the gap parameter on the chemical potential in 
the mean-field approximation. These results are analytically expressed in terms of elliptic integrals 
of the first and second kind. We discuss the form of these coefficients in the extreme BCS and 
BEC regimes and around the unitary limit, and compare with recent work by other authors. 



I. INTRODUCTION 



The last fewyears have witnessed a renewed interest in the physics of the BCS-BEC 
crossover [H, 0, y], partly motivated by the availability of tunable interactions in the realm 
of interacting Fermi gases [4|]. Recent experimental work has shown evidence 

for condensation of fermionic atom pairs, suggesting the formation of a fermionic superfluid. 
From the theoretical point of view, the qualitative description of the BCS-BEC crossover 
has been based on the mean-field theory of Leggett [2J and its extension to finite tempera- 
ture by Nozieres and Schmitt-Rink [Io| and Sa de Melo, Randeria and Engelbrecht [ll[. In 
this description of the superconducting system, the effective action in terms of a complex 
order parameter which couples to the pairing field plays a central role. Some recent devel- 



opments [12|, [13|, [14|, [15j in the crossover problem beyond mean field theory have improved 



our understanding of the equilibrium state at a quantitative level. In particular, Diener et 
al. ill, by computing the complete quadratic part of the effective action, have obtained 



the correction to the mean field result which arises from the integration of the Gaussian 
fluctuations, finding excellent agreement with calculations based on quantum Monte Carlo 
techniques in the unitary limit where the scattering length a — > 00. 

The effective action can also used, in principle, to derive an effective Lagrangian which 
captures the low-energy behavior of the system in terms of the Goldstone mode phase of 
the order parameter. At zero temperature, one expects that an expansion of the effective 
Lagrangian in derivatives of the Goldstone field can be used in order to study the low-energy 
behavior of the system. The lea ding order (LO) in this expansion was evaluated by Greiter 
et al. [161] and by Aitchison et al. [171 some years ago, and since then, there have been various 



microscopic derivations (l8l . 19 . 20 . 21, 22| of effective models, some of which (l9l . 20 ] have 
included more or less explicitly some derivative corrections. 

In a recent work, Son and Wingate 23f| have systematically studied the form of the ef- 
fective Lagrangian in the unitary limit at the next-to-leading order (NLO) in the derivative 
expansion, when the effective theory is formulated in terms of the Goldstone mode coupled 
to external gauge and gravitational fields. At unitarity, it turns out that, besides general 
coordinate and gauge invariance, the theory exhibits conformal invariance [241], which puts 
constraints on the form of the NLO Lagrangian and restricts to two the number of indepen- 
dent NLO parameters. 

In this paper, we extend these studies. We evaluate the effective parameters of the 
NLO Lagrangian in the unitary limit at the mean-field level, and also obtain the most 
general form of this Lagrangian away from this limit, where the symmetry under conformal 
transformations is not realized. By computing all the necessary functions at the one-loop 
level in terms of elliptic integrals, we have obtained the simplest approximation to the low 
energy effective theory for the whole crossover region at zero temperature. 

We find that, at unitarity, the effective Lagrangian is specified by two constants, but its 
actual form differs from that given in Ref. 23||, and includes a new contribution (didjO) 2 
of higher spatial derivatives of the Goldstone mode, while the NLO contribution of the 
external trapping potential proportional to V 2 V ext is absent. We show that these features 
are a necessary consequence of the conformal invariance of the NLO field equations. As an 
application, we derive the energy density functional in the unitary limit, and compare it 
with the computation of Rupak and Schaffer [251 ]. which is based on an epsilon expansion 
around d = 4 — e spatial dimensions and is, to our knowledge, the only one in the literature. 
Altough the coefficients computed by these two methods show discrepancies of the order of 



2 



30% in general, we find a surprisingly good agreement for the coefficient of the quantum 
pressure. 

For the whole crossover region, we obtain the NLO Lagrangian in terms of four functions 
which are given in closed form in terms of elliptic integrals. The BCS and BEC limits as 
well as the near unitarity limits of the NLO Lagrangian are worked out in detail. In the 
BEC limit, we recover the known features of the hydrodynamic description of superfluidity 
at zero temperature. 

The plan of this paper is as follows. In Section [TT] we formulate the problem in the 
framework of linear response theory and derive a linearized equation for the Goldstone field 
in terms of derivatives of response functions. In Section II I II we show how to construct a 
Lagrangian, including second order time derivatives of 9, by considering all the available 
Galilean invariants consistent with required general properties. Then we present a careful 
procedure of reduction, and show how to use the LO field equation in order to eliminate 
undesired higher order time derivatives without changing its perturbative contents. We also 
argue the need to compute two three-point functions in order to determine all the coefficient 
functions in the effective Lagrangian. In Section [IV] we present an analytical expression for 
the thermodynamical potential at the one-loop level, and the analytical expressions of the 
NLO coefficients in the two (BCS and BEC) limits and near unitarity. In Section [V] we 
compute the energy density functional in the unitary limit, and compare our result with 
other approaches. Section [VH gives our conclusions. Details of the calculations as well as 
additional material on the invariance properties under conformal transformations in the 
unitary limit are given in a series of four appendices. 



II. DERIVATIVE EXPANSION OF THE RESPONSE FUNCTIONS AT T = 

The system is conveniently described by the BCS Lagrangian density in terms of a Nambu 
spinor field and a complex field which decouples the short-range interaction in the Cooper 
channel 



id t + r 3 -^V 2 + Ut 1 + it 2 )A + Ut 1 - it 2 ) A* 



A*A, (2.1) 

9A 



where ^ = (ipt, ip\) and r* is the corresponding Pauli matrix. The bare coupling parameter 
g\ depends on the details of the short range interaction and on a regulator that truncates 
some loop integrals at some very large scale A. As usual in the framework of effective field 
theories, the result of adding the appropriate loop corrections to this bare coupling will be 
matched with the measured low energy scattering properties encoded in the s-wave scatter- 
ing length a. Other physical parameters such as the s-wave effective range or the p-wave 
scattering length are related to operators with more derivatives in the effective Lagrangian, 
and hence, they are generically subleading in the expansion in powers of Rd, where R is 
some length setting the size of the interaction region. The Lagrangian is invariant under the 
U(l) symmetry of phase independent spacetime-independent transformations ^ — > e lT 
A — > e l2e A, which is spontaneously broken down to H = Z 2 below a critical temperature. 
In order to compute the effective action for the resulting gapless collective mode it is conve- 
nient 26[ to express the fields as a U{1) transformation acting on fields \l/ and a which do 



not contain the Goldstone mode 

q?( x ) = e irH ^{x), A(x) = e l29( - x \A + a(x)), (2.2) 



3 



where the constant amplitude A and its fluctuation a(x) are real numbers. Since we shall 
hereafter use the fermion fields the tildes will be dropped to simplify the notation. With 
this choice, the Lagrangian involving the couplings between the fermion field and the other 
fields becomes, after an integration by parts, 

Cj = -$t r 3^ ( dt Q + _l_(v#) 2 + V ext J - — (^ f W - Vtf f tf) • V9 + #V*<7, (2.3) 
\ 2m J 2mi 

where V ext is an external potential. The quadratic part £ 2 of the Lagrangian including the 
chemical potential coupled to conserved particle density \|/V 3 \I/ is given by Eq. (12.11) with 
the replacements V 2 /(2m) — > V 2 /(2m) + //, A — > A and A* A — > a 2 . 
Under an infinitesimal Galilean transformation 



x 



x 



x + vt, 



t^t' = t, 



(2.4) 



the Goldstone field changes inhomogeneously as 



56{t,x) 



-vt ■ V6 1 + mv ■ x. 



(2.5) 



Once the auxiliary field a has been eliminated, the invariance of the effective action for 9 
requires a dependence through the scalar quantity 



X = n - d t 9 



2m 



ext 5 



5X = -vt ■ VI, 



(2.6) 



and some appropriate derivatives of X and V0. According to the power counting scheme of 
Ref. [23(| 1 . these quantities are 0(p°). This can also be seen by noting that, at equilibrium, 
the gradient of the phase v s = m~ 1 'V9 is the constant superfluid velocity, and d t 9 goes 
with the chemical potential. To next-to- leading order, we have the following 0(p) galilean 
invariant derivatives 

V9 



V z 9, 



d t + — -V 

m 



X, VI, didjO, 



(2.7) 



where we have only written derivatives of the quantities X and V0, which are coupled to 
the particle density and the current. The signature under time reversal (t — > —t, 9 — > —9) 
of the two scalar terms (and the tensorial term) is —1, and hence there is no contribution 
of 0(p) to the effective Lagrangian. In the next Section we will list all possible scalar terms 
of 0(p 2 ) potentially contributing to the effective Lagrangian. 

A possible way to derive the effective Lagrangian for the Goldstone mode is to compute 
the appropriate linear response at low frequency and momentum for an external perturbation 
given by a Hamiltonian 



tfVtf ta0 + Kxt) + ^(^ f V^ - W^) • W - tfVtftr 

2mi 



H ext (t) = J d 

= J d 3 x [n(d t 9 + 1/ C xt) + 3 P ■ V# - Vtftr] 



(2.8) 



A quantity containing N[dt) time derivatives, N[di] spatial derivatives and N[9] powers of the Goldstone 
field is counted as 0{p N ), where N = N[d t ] + N[di\ - N{8}. 
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where j p is the current density operator in the absence of V# and the total current operator 
is 3 — 3 V + ^T 3 ^m~ 1 'V9. The induced changes to be computed are 6{n(Q)), S(j(Q)) 
and the change in the expectation value of the pairing field <5(\I/V 1 \I/) in terms of 9{Q) 
and cr(Q), where Q = (q,ou). By combining the conservation of the particle number with 
the gap equation for 5(\E'V 1 \[ / ), one can obtain a single linear equation of motion for 9 
and, consequently, a quadratic Lagrangian to be consistently matched with the form of the 
Galilean invariants listed above. 

In what follows, we use Xab{Q) for the Fourier transform of the retarded response function 
Xab(X,X') = -i([A(X),B(X')])9(t - f). For the Hamiltonian given in Eq. ([ZED, the 
required set of linear response equations is 

S(n(Q)) = Xnn(Q) (-iu9(Q) + V cxt (Q)) + Xn %Q)^(Q) ~ Xm(QMQ), (2-9) 

5(j k (Q)) = x%{Q) Hco9(Q) + Kxt(Q)) + xfj(QWo(Q) - x k n (QMQ), ( 2 - 10 ) 
a^V 1 *) = xm(Q) (-ico9(Q) + v cxt (Q)) + Xi-mq k o(Q) - xu(.QHQ), (2.11) 

where there is a summation over upper repeated indices. It is useful to collect some symmetry 
properties that follow from the behavior under time reversal and parity: 

X%{Q) = Xn){Q) = Xjn(Q)q\ (2.12) 
Xn(Q) = Xi*(Q) = Xji(Q)q k , (2.13) 

Xnl{Q) = Xln(Q). (2.14) 

When the operators A and B have the same (opposite) signature under time reversal, 
Ihixab(Q) is odd (even) in uj. This implies that Re Xjn(Q), Re XjiiQ) are °dd in uj, while 
R e Xni(Q), ReXnn(Q) and Rexn(Q) are even functions of uj. The current-current response 
function can be written as 

xfj(Q) = XL(Q)q k q l + Xt(Q) {S kl - q k q l ) , (2.15) 

where Rcxl,t(Q) are even in uj. At q = 0, they must obey the usual sum rules giving the 
static responses at superfluidity 

XL (0) = lira ^ ImXL( ^ } = M, (2.16) 

XT (0) = lim f - ImXT(g ^ } = {n) ~ U \ (2.17) 

where n s is the superfluid particle density. 

The results of this Section do not depend on the specific approximation used to compute 
the response functions, but in order to gain some physical insight on them, we write the 
kind of integrals to be computed at the one-loop level 



(2.18) 



J W { (2k + q) k {2k + q) 1 ) 
Here r° is the identity matrix and Q is the Nambu-Gorkov Green's function 

g <^> ^ + ^J A ° ' (2 ' 19) 
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with = €fc— = k 2 /2m—fi and = v/£|~+ Ag. Note that, together with the contribution 
coming from [jf p ,j p ], the current-current response function X^-(Q) includes the longitudinal 



piece (n)m -1 ff, where (n) is the total particle density. The integration over the imaginary 
frequency produces the denominators u + is ± (i?fe + Ek +q ), which for small frequency 
\u\ < 2Ao do not contribute to the imaginary part of the response function. In addition, the 
expansions about u = and q = are well behaved, in marked contrast to the expansions 
of the denominators u + ie± (E^ — Ek +q ) which arise from Landau damping at T ^ 0. Thus, 
at T = the needed response functions are regular real functions near Q = 0. After some 
calculation, the one-loop expressions for Xab(q = 0,u = 0) become 



Xm(0) = - / — ^ (2-21) 



d 3 k A q 
(27r)3 El 

X[V(0) = - r^lr, (2.22) 



(2tt)3 i?3 

4(0) = (2.23) 

and we see that static current transverse response Xt(0) vanishes which, according to 
Eq. (12.17j) . shows that the entire system is superfluid, (n) = n s , at T = 0. 

The linear ultraviolet divergence in xn nas been regulated by the cut-off A. As it is 
well known, the renormalization of this divergence is performed by the substitution of the 
bare coupling constant g\ in terms of the s-wave scattering length a. Putting together the 
leading order interaction of the form g^ 1 and the first correction in the vacuum, one obtains 
the measured coupling constant g = —Anam' 1 , 



1 1 f A d 3 k 1 



9 9a J (27r)3 2e fc : 



^2.24) 



and this relation allows the determination of g\ in terms of a. As all the calculations will 
be A independent, it is convenient to define the 'renormalized' Xn as 



As mentioned above, at T = all response functions are real in the region of small Q and 
their lowest-order derivatives at Q = must determine the form of the effective Lagrangian 
at next-to-leading order in the derivative expansion. We next show how to compute this 
Lagrangian. From the gap equation 

S^t 1 *) = —cr{Q), (2.26) 
9a 

and the equation for the response S^W 1 ^), we can write the change <r(Q) in the local 
amplitude in terms of 6(Q) and V ext (Q). If we replace the result for a(Q) in the continuity 
equation 

-u5(n(Q)) + q k S(j k (Q)) = 0, (2.27) 
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and expand in powers of uj and q, we find the following equation for 9(Q) 
= (hu 2 + b 2 q 2 + b 3 u 2 q 2 + b 4 u 4 + b 5 q 4 + . . .) 9{Q) 

+ (61 + + huj 2 + ...) iuV ext (Q), (2.28) 



where 



62 = -Xl(0) = -M, (2.29b) 

m 

6 3 - hit" - ~ o , — -*r: hit _ /o , — ^ ^ ' ^-2ycj 



9c; <9g 2 9a; 2 2 + gxn V da; 9q 2 J (2 + gxn) 2 <9<? 2 

dXnn . ZgXnl ^XnY S^Xnl ®X 



7 "AMI _ ~^A.rtl "A. Til _ ^ A.nl /<-, nn n 

4 ■ &0 2 + (2 + g Xll ) 2 du;^ 1 ' J 

h = (2.29e) 

, <9X> _ ^frm _ ggXnl ( 1 ^Xjl _ gXrU \ _ ^Xnl gXll /« 2qf x 

6 " dg 2 2 + gxu \2 du ' «9g 2 J (2 + g X u) 2 dq 2 ' 1 ' J 

All these parameters are computable from the appropriate two-point retarded functions and 
their derivatives. As the gap equation in terms of the thermodynamic potential f2(/i, A ) at 
T = 

— X A = 2.30 

implicitly determines the function A (/x), the coefficients depend on the chemical po- 
tential and, parametrically, on the scattering length a. 



III. THE HIGHER-ORDER AND THE REDUCED EFFECTIVE LAGRANGIANS 
A. The higher-order effective Lagrangian 

We now write an effective Lagrangian for the Goldstone mode which leads to the above 
equation of motion for 9. This Lagrangian contains higher-order time derivatives, which 
will be dealt with in the second part of this Section. Once the auxiliary field a has been 
eliminated through the use of the gap equation, the NLO Lagrangian is a linear combination 
of all the independent scalar operators of 0(p 2 ) constructed from derivatives of X and V0, 
with coefficients which are functions of X, the only scalar of 0(p°). Note that, as mentioned 
above, 0{p) scalar terms such as ri(X)V 2 # s and r 2 (X)X( are excluded by invariance under 
time reversal. 

There are only five independent NLO terms, (VX) 2 , (V 2 9) 2 , {didjOf, X 2 , and X t V 2 9, 
where 

X t = d t X + to -1 V0 • VX . (3.1) 

Other candidate terms can be shown to be dependent on these five. For instante, we have 
excluded a term /i(X)V 2 X , which becomes proportional to (VX) 2 after an integration by 
parts. Another potential contribution of 0(p 2 ), (d t + m^ 1 'V9 ■ V)V 2 #, which by 

(d t + m~ 1 V9 ■ V)V 2 9 = -V 2 X - VV cxt - m" 1 ^^) 2 
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is equivalent to V V ext , must also be excluded, since the identity 

G\X)X t W 2 6 + G'(1)(VX) 2 - G(X)VV ext - -j-G(X) [{d^ef - (V 2 6) 2 ] = 
d t {G(X)V 2 6) + V • G(X)V 2 6 + G(X)Vlj (3.2) 

shows that this term is in fact redundant. Thus, the most general Lagrangian up to next- 
to-leading order in the derivative expansion is given in terms of six coefficient functions 

C = P{X)+g l {X){VXf + g 2 {X){V 2 ef + g,{X){d i d j e) 2 + g,{X)X 2 t +g b {X)X t V 2 e i (3.3) 

to be evaluated from the low-energy behavior of the response functions. 
The linearized field equation following from this Lagrangian has the form 

+2 [g 2 (i,) + gM) + [29,0m) - gM) + 2g^) + (3.4) 

Comparison with Eq. (12.281) yields the relations 

9i = -y + h, (3.5a) 
92 + 93 = y, (3.5b) 

94 = y, (3.5c) 

95 = ~h + b 6 , (3.5d) 

and 

P"(/i) = b u (3.6a) 

P'^) = -mb 2 = (n). (3.6b) 

These expressions match the coefficients of the effective Lagrangian with the response 
functions. Note that from the response functions we can only determine the sum 
g 2 (X) + g 3 (X), but not the individual functions. This is not surprising since, in the 
quadratic approximation to the Lagrangian, the terms (V 2 #) 2 and (didjd) 2 are not inde- 
pendent, but differ by a total derivative. Hence, in order to evaluate g 2 and g 3 separately, 
one must resort to the computation of three-point functions. Details of this computation 
at the one-loop level are given in Appendix [Dj where we find the remarkably simple result 
g 3 (X) / g 2 (X) = 2. This implies: 

g 2 (X) = h 5 , g 3 (X) = \b 5 . (3.7) 
6 3 

Regarding the expressions for the leading order coefficients bi(fx) and b 2 (fj,), it is possible 
to check their consistency without actually computing them. As the authors of Refs. [l6, 23| 
have shown, the leading order effective Lagrangian at T = is a function of the invariant 



8 



X which precisely coincides with the pressure P(/x) as a function of the chemical potential. 
Obviously, the expression for 6 2 i n (12.29bj) satisfies this requirement. To see the agreement 
of b\ in f)2.29ap with P"(/i), we note that the pressure is given by P{p) = — A (//)), 
where Ao(/x) satisfies the gap equation ( 12.301) . The thermodynamic potential can be written 
as the sum of a 'tree' level renormalized piece and the full quantum contribution $ from the 
perturbative expansion 

!i ("' A ») = | + < t(ft A ). (3.8) 



V 



By using A' (fi) = - (d 2 Vt/dA 2 ) d 2 n/dfidA , a direct evaluation of P"(/x) yields 



d 2 Q f d 2 Q 



<9 2 f] 

dAj 



-1 



(3.9) 



=a (m) 



which agrees with the expression for b\ when one identifies properly the static response 
functions with the susceptibilities 



Xnn(0) 
Xnl(0) 



V 



V 



xn(o) = v- 



dfi 2 ' 

dfidA ' 
_ l( 9 2 ft 2 



9A 2 



(3.10) 
(3.11) 
(3.12) 



Note that the condition of thermodynamic stability implies that b\ must be positive. This 
analysis shows the important role played by the response functions of the pairing field 
vl/V 1 ^ in the construction of an effective Lagrangian satisfying known general properties. 
Some recent work 22] has overlooked this point. 



B. Reduction of the higher-order Lagrangian 

To proceed further, we must find a "reduced effective Lagrangian" without higher-order 
time derivatives, but perturbatively equivalent to Eq. (13.31) . The LO field equations can be 
used, in principle, to eliminate terms with higher order time-derivatives in favour of terms 
with spatial derivatives but, in doing so, one must be very careful that the perturbative 
content of the original Lagrangian is preserved. In this regard, it is important to note that 
the LO equations 

x ' + ^§) v2e = ^ < 3 - 13 > 

are satisfied only up to terms of 0(p 3 ) coming from the NLO Lagrangian. 

The term proportional to X 2 , which gives rise to fourth-order time derivatives of 9 in the 
field equation, can be eliminated by adding and subtracting g±(X) times the square of the 
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LO field equations. This yields 



C = P(X) + g 1 (X)(VX) 2 + g 3 (X)(d i d,ef 



+ 



+ 



g 2 (X)-g 4 (X) 
g 5 (X)-2g±{X) 



P'{X) 
mP"{X) 

mP"{X) 



(V 2 



x t v 2 e + 9i (x) 



x t 



mP"(X) 



l 2 



v 2 e 



(3.14) 



Now, the term proportional to the square of the LO field equation, when evaluated on a 
perturbative solution, is of 0(p 6 ), and thus highly suppressed. Furthermore, one can easily 
check that its contribution to the field equation is of 0(p 5 ). It can thus be safely dropped 
and we are left with the following reduced Lagrangian 

I = P(X) + 9l (X)(VX) 2 + g 2 (X)(V 2 6) 2 + g 3 (X)(d t d,6) 2 + g 5 (X)X t V 2 6. (3.15) 

We still have to get rid of the second-order time derivative in the last term of £nlo But 
in this case we can not use the leading field equation 



g h {X)X t V 2 9 = g 5 {X) X t + 



mP"(X) 



v 2 e)v 2 e-ux) 



mP"(X) 



(V 2 



(3.16) 



to this end since, even though the numerical correction to the classsical effective action 2 
introduced by dropping the first term in the RHS of Eq. (I3.16P would be of 0(p A ), and 
thus acceptable, the corresponding correction to the field equation would be 0(p 3 ), which 
is of the same order as the NLO contributions to the field equations. This can be checked 
explicitly by computing the linearized correction to the field equation coming from the first 
term of the RHS of Eq. fl3~TBl 



- 2£ 5 (/i)V 2 



d 2 e dv„ 



dt' 



+ 



dt mP"(fi) 



dt 



(3.17) 



While, by Eq. f l 3 . X 3 [) . the first term is 0(p 5 ) on a classical solution and can be dropped, the 
second term yields a contribution of 0(p 3 ) which produces an unacceptable change in the 
perturbative field equations. 

Fortunately, we can instead use integration by parts, which is guaranteed to exactly 
preserve the numerical value of £nlo and the perturbative field equations. With the help 
of identity (13.21) . we perform the replacement 

Mx)x t v 2 e - -g,{x){vx) 2 + G 5 (x)vV ext + -g 5 (x) [{d^ef - {v 2 ef] , (3.18) 

where G' 5 (X) = g 5 (X). Thus, the reduced effective Lagrangian becomes 



c = p(x) + h{x){vxy + f 2 (x)(v 2 ey + h{x)(d t d 3 ey + u(x)v 2 v e 



ext ; 



(3.19) 



By this we mean the effective action J d 3 xdtC evaluated on a classical solution. 
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where the coefficient functions are given by 



f x {X) = 9l (X)-G' 5 (X), 



(3.20a) 




2 



1 = 



G 5 (X) 



(3.20b) 



m 



f 3 (X) = 93 (X) + -G 5 (X), 
m 

U(X) = G 5 (X), 



(3.20d) 



(3.20c) 



with 



G' 5 (X) = g 5 (X)-2 9i (X) 



mP"(X) ' 



(3.21) 



The constant of integration in G^(X) is irrelevant, since it enters the effective Lagrangian 
as the coefficient of a total divergence, as follows from Eq. (13. 2p when G(X) is a constant. 

Note that the term proportional (didj8) 2 , which has not been considered previously in the 
literature, is an essential ingredient of the effective Lagrangian and can not be eliminated 
without introducing unacceptable changes in the NLO field equations: Although we could 
use Eq. (13.21) to eliminate this term in favour of the other three invariants, in doing so we 
would reintroduce the term proportional to X t 'V 2 9 which, as we have seen, should not be 
eliminated through the use of the LO field equations. 

An alternative way to understand this issue is by noting that using the LO field equation 
in Eq. (13.161) would be equivalent, up to 0(p 4 ) terms arising from the second variation of 
P(X), to performing the field redefinition 3 9 = 9 + 59 



The variation of P(X) under this redefinition cancels the last term in Eq. (13.151) and we are 
left with an effective NLO Lagrangian depending on only three coefficient functions. Since 
a field redefinition can change neither the value of the action evaluated at its extrema (the 
classical action) nor the result of any functional integration (the generating functional), it 
is clear that, for some applications, one can use an efffective NLO Lagrangian with only 
three coefficient functions. However, the field redefinition (13.221) involves a change of 0(p 2 ), 
making 9 and 9 non-equivalent at NLO. Thus, with hydrodynamical applications in view 
where the field 9 has a concrete physical meaning — through the relation to the superfluid 
velocity v s = m~ 1 'V9 — one is forced to keep all four NLO coefficient functions in (13.191) . 
Only this way can we preserve, at NLO, the interpretation of 9 as (1/2) the phase of the 
condensate, the canonical structure of the theory and the perturbative field equations. Note 
that the situation is different with the use of the LO field equation to eliminate the term 
proportinal to X 2 in Eq. (I3.14p . which is equivalent to a field redefinition with 



In this case the difference between 9 and 9 is of 0(p 4 ) and they are equivalent at NLO. 



3 We are indebted to the authors of Ref. [23| for this observation. 
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P"(X) 



V 



(3.22) 




(3.23) 
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IV. EFFECTIVE PARAMETERS IN THE MEAN-FIELD APPROXIMATION 



The effective theory of the previous Section contains a set of functions depending on 
the chemical potential and the scattering length once Ao(/i) is inserted. Calculating these 
functions in the mean field approximation is relatively straightforward but tedious and, as 
explained in Appendix A, all the integrals ( 13. 15ft are computable in closed form in terms of 
complete elliptic integrals of the first and second kinds; the results are given in this Section 
and in Appendix O Other physical quantities, such as the length for pair correlation, 
have been expressed in terms of elliptic integrals in Ref. 27|, but here we will focus on the 
thermodynamical potential and the coefficients fi(X). Henceforth the particle density is 
written as (n) = k F /37T 2 , where kp is the Fermi wave vector. 



A. Leading order results 

Here, we briefly collect the results for the pressure and the static, zero-momentum sus- 
ceptibilities in the mean field approximation. We start with the one-loop thermodynamic 
potential [l5[ 



n 
v 



Ana 



d 6 k 



1 Aq 

-Sfc — Cfc + ^ — - — 



(4.1) 



which, as mentioned above, can be differentiated to yield Eqs. (12.201) . (12.211) and (12.221) in 
the same approximation. The integral can be done anlytically, giving 



n 
v 



mAl 4m 3/2 |/i| 5/2 



Ana 



157T 5 



3a' 



-si 



ign(/i) + v / l + « 2 E(-7) 



sign(//)(l + a 2 ) + (3a 2 - 1)^1 + a 2 



K(-j) 



-sign(/i) + vT+a 2 



(4.2) 



where K(n), E(n) are the complete elliptic integrals of the first and second kind respectively 
(see Appendix A for details and notation), and a 2 = A 2 ,/^ 2 . The parameter 7 is 



7 = 

The susceptibilities are given by 

Xnn(0) 



+ a 2 + sign(/i)] ' 



a' 



Xnl(0) 

Xn(0) 



m 3/2|^|l/2 



nr 



3/2 A 



sign(/i) + Vl + a 2 [-£(-7) + if (-7)] , 
1 



7T 2 |/i 



1/2 



m 3/2|^|l/2 



-sign(/z) + Vl + a 2 
1 



7T 



a 



■VT+ 



a" 



— sign(/i) + Vl + a 2 3a 2 Vl + a 2 E(--f) 



2sign(/i)(l + a 2 ) + (2 + 3a 2 )Vl + a 2 K(-j) 



(4.3) 

(4.4) 
(4.5) 



(4.6) 
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The expression for the thermodynamic potential (14. 2p . the gap Eq. (I2.30p and the number 
equation 

{n) = -van (4 " 7) 

summarize the properties of the crossover at T = in the one-loop approximation. From 
these one can find the quantities A /ep and \ijep as functions of the parameter (kpa)' 1 , 
where ep = k F /2m is the Fermi energy. The energy density of the ground state is given by 
e = fi(n) + V- l tt. 

Next we collect results, some of then well known, which follow easily from Eq. (14.21) for 
specific ranges of the parameter (kpa) -1 . 



1. Near unitarity 

Here we present the results for the gap parameter, the chemical potential and the ground 
state energy per particle near unitarity, expressed as power series in (/c^a) -1 : 

A 0.6368 0.0959 , , 

= 0.6864 + — + _— + ..., 4.8 

ep kpa {kpa) z 

a 0.7401 0.5150 , 

^ = 0.5906 -— 77r^ + ---> 4-9 

ep kpa {kpa) z 

e 3e F / ncnn „ 0.9251 0.8582 \ /A 1 . 

- r -(0. 590e - — - — + ...j. (4.10) 

The numerical coefficients have been obtained by simultaneous power series solution of the 
gap and number equations using the analytic solution (14. 2p for Q. It is also easy to write a 
few terms for the pressure at large scattering length 

2 1/2 3/2 

PU) = 0.0842 m 3 V /2 + 0.1075^- + 0.1274 m f + O.lOOe^ + • • • (4.11) 

a a 1 cr 

2. BCS limit 

In the BCS limit (kpa)' 1 — > — oo. By using the following asymptotic expansions for the 
complete elliptic integrals, valid when z — > — oo 

~ ^P + ^^ + OK-^M-.)), (4-12) 
- (-») ^ + + 3 } + ° «-)- 5/2 M-»)) ■ (4-13) 

we obtain the solution to the gap equation a = 8e~ 2 exp(7r / 2y/2mfia)) . This produces the 
pressure 

2 5/2 3/2 5/2 

P(fj) = 9 (1 + 60 exp(-4 + vra- 1 (2m/i)- 1/2 ) + ...), (4.14) 

157T Z 
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and 



%-Ke- i e K/kFa 
1 + 1 + 



1 + 



A 2 
44 



2 + In 



8e F 



(4.15) 
(4.16) 



We find that the ground state energy per particle and the fermion density are given by 



e 

V) 

(n) 



IT 

3^ 



3Af 

8e F 



+ ... 



1 + 



164 



1 -2 In 



8e F 



(4.17) 
(4.18) 



3. BEC limit 

The other extreme regime, the BEC limit (kpa)^ 1 — > 00, is obtained when \x is close to 
— l/2ma 2 . In this regime |a| is small and we can use the Maclaurin series for the elliptic 
integrals. The thermodynamic potential is thus 



and one finds 



n 

V 



e 



mA 2 m 3/2 ( _^l/2 A 2 



III 



3 / 2 A^ 



4na 



2ma 2 



8ji = e F 



647r v / 2(-/i) 3 / 2 



+ 



2k F a 



+ 



{^k F a) l l 2 



+ 



+ 



kpCL 



(k F a) 2 3tt 

The solution of the gap equation yields a power series in ma 2 5fi for Ao 

5 

H — ma 2 8jji + 0((ma 2 5n) 2 ) 



2v% 



a\/m 



which can be used to write the pressure as 



P 



m5/i2 +^^ + ... 



2-na 



An 



(4.19) 

(4.20) 
(4.21) 
(4.22) 



(4.23) 



(4.24) 



Eliminating the chemical potential leads to a pressure P(n) = (2m)~ 1 TTan 2 + . . .. Noting 
that iris = 2m and tib = n/2, a comparison with the pressure of a weakly interacting Bose 
gas in the lowest-order approximation P(ub) = 2%aBm B /1 n B , yields the mean-field result 
for the scattering length between bosons, as = 2a. But, as Diener, Sensarma and Randeria 
have recently shown 15j , these mean-field results in the BEC limit are poor approximations 
to the results which are obtained when the contribution from the quantum fluctuations is 
included in the computation. 
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B. Next-to-leading order results in the one-loop approximation 

The derivatives of the one- loop response functions at Q = are collected in Appendix O 
From these expressions we can determine the coefficient functions fi(X) in the effective NLO 
Lagrangian 

£nlo = fi(X)(VX) 2 + f 2 (X)(V 2 6) 2 + hWididjO) 2 + h(X)V 2 V ext . (4.25) 



1. Unitarity 

It is possible to express the complete elliptic integrals in the unitary limit through two 
useful formulas obtained by simultaneous solution of the pressure P = com 3 / 2 /z 5//2 and gap 
equations. These are given by 



i5ttV-i + yrr 



a 



2 



K (~^ = TrTZ 2^ c o> ( 4 - 26 ) 

4(1 + cr) 

E{-l) = ■ en. (4.27) 
4Vl + a 2 V-l + VI + a 2 

Substituting these expressions into the derivatives of Appendix C produces the remarkable 
result that G' 5 (fi)—0. This is due to the fact that, when the coupling g — > oo, the expression 
for this quantity is proportional to the gap equation, which gives rise to the exact cancellation 
of G' 5 (fi). The complete list of coefficient functions that we find is 

„ . , 35c m x l 2 m 1 / 2 

fM = -^^72" - -0-0153-^, (4.28a) 
hip) = -J^T^^m - -0-0035^—, (4.28b) 

fM = 0, (4.28d) 

where the numerical values are obtained by substitution of the one-loop numerical values 
a(n) ~ 1.1622 and c ~ 0.0842. Noting that f 3 /f 2 = -3, our result for the NLO effective 
Lagrangian at unitarity can be written as 

C NLO = Cl m l l 2 X~ l / 2 {VX) 2 + c 2 m- l l 2 X 1 ' 2 [{V 2 6) 2 - 3(^5^) 2 ] , (4.29) 

where C\ ~ —0.0153 and c 2 ~ —0.0035. This is one of the main results in this paper. 

One might wonder whether the cancellation of the f± coefficient and the ratio jzj f 2 — — 3 
found here are mere accidents of the mean field approximation. Actually, these are exact 
consequences of the conformal invariance displayed by the system at the unitarity limit. In 
other words, any approximation scheme that respects conformal invariance must necessarily 
yield a result of the form given by Eq. (14.291) . This is explained in Appendix [Bj where the 
reason for the discrepancy with the form of the effective action given in [23[ is also analyzed. 
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2. Near unitarity 



The following expansions are valid for the NLO coefficient functions near unitarity 



/ 2 (a0 



0.0153m 1 / 2 0.0014 0.0016 

a fx a 2 m}/ 2 fx 3 / 2 



1/2 



1/2 



0.0035/i 

m l/2 

0.0004 



a 2 m 3/2^1/2 



+ 



+ 



0.0023 + 0.0022 In — 
0.0022 



+ 



1 

am 



a 3 m 2 fx 



+ • 



1/2 



0.0104/i 

m l/2 

0.0097 

a 2 m 3/2^1/2 

0.0022 , / u 
In — 

a \/i 



0.0113 + 0.0022 In 



JL 

/'o 



1 

am 



0.0094 

a 3 m 2 /i 



+ • 



0.0037 



0.0004 
+ — + ..-, 



a 2 m i/2^i/2 a 3 m\x 
where /xq is an arbitrary scale (ln/xo is multiplied by a total divergence). 



(4.30a) 



(4.30b) 



(4.30c) 
(4.30d) 



3. BCS limit 

The determination of the leading behavior of the coefficient functions in the BCS limit is 
more involved. By using the asymptotic expansions for the complete elliptic integrals and 
the perturbative solution of the gap equation a(fx) ~ 8e~ 2 exp(7r/2a v / 2myu), one finds the 
lowest-order approximation 



0Xnl(O) 



0(a 3 (lna) 



(4.31) 



2 + ^Xn(0) 2 

to be used in the coefficients of some derivatives in the asymptotic expressions for the bi 
coefficients. The pressure (14.141) gives 



2/x ix 2 / a(/i) 
— — a(fxy ha 1 



2 + In 



a(fx) 



mP"(/x) 3m 6m 
and the substitution of these results into Eqs. (I3.20p produces 

1 



36V2tt 2 n 



1/2 



and 



/i(A0 
/a 00 
/s(/*) 

AO*) 



(hiJa^V^mV 2 
24v^tt 2 /U 1 / 2 

2 3/2 



+ . . . 



288^2 mV2aV /2 
1 



+ 0(a 4 (lna) 4 ), 



+ ..., 



1/2 



1357r 2 a;(^) 2 m 1 / 2 



2 3/2 



1/2 



457r 2 a( / u) 2 m 1 / 2 
1 



+ . . . = 

+ 



192 v / 2m 1 /2 a 2 /U 3/2 
exp(4-7ra- 1 (2m / u)- 1 / 2 ) li 1 / 2 

m l/2 
+ • • 



2160\/27r 2 
exp(4-7ra- 1 (2m/i)- 1 / 2 ) ^/ 2 

V2" 



+ 



720v^ 



7T 



m 



144V2m 1 /2 a 2 At i/2 



(4.32) 

(4.33) 

(4.34a) 
(4.34b) 
(4.34c) 
(4.34d) 
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It is worth noting that f 4/ {mf = 0(a 2 (ln a) 2 ), which shows that can be safely ne- 
glected. Hence, an important feature of this regime is the dominance of the fi and terms 
in the next-to-leading order effective Lagrangian. The coefficients 64 and 65 are at the origin 
of the leading behavior of ^2,3, while 63 and 64 govern the expression for j\. 

Note that these expressions can not be trusted for /c^a arbitrarily small, because then 
Ao — > 0, and the condition for the validity of the derivative expansion \u)\ < 2Ao cannot be 
satisfied. Indeed, when a —>■ 0~ the Goldstone mode can hardly be considered a propagating 
mode due the arbitrarily small two-particle states threshold, and this renders this effective 
field description meaningless. 



4. BEC limit 

The BEC limit can be obtained by substituting the gap parameter (14.231) into the appro- 
priate derivatives of the response functions and then expanding in powers of 5fi. We obtain 
the following leading behavior for the coefficient functions 
~ 777/7 

G'M = -— + 0(8»), (4.35a) 

4o7T 

I th n 

= + 0M = - iMi + Wri + (435b) 

fM = -7^— + 0(^ 2 ), (4.35c) 
9o7rma 

f ^-mL + ^- < 435d ) 

UP) = ~ + 0(S^). (4.35e) 

4o7T 

These results reveal that, in this limit, G' 5 (n) / f\(fi) = 0(ma 2 5fi), and the leading derivative 
contribution comes from the f\ term of the effective Lagrangian. Now, the coefficient 63 is 
at the origin of the leading behavior of f±. Note that the constant coefficients and ^3 play 
no role in this limit, because an integration by parts produces their mutual cancellation. 



V. DENSITY FUNCTIONAL IN THE UNITARY LIMIT. RELATION TO OTHER 
APPROACHES 

From the effective Lagrangian one can easily derive the energy density £ depending on 
the number density and spatial derivatives of the Goldstone mode. The corresponding first- 
order equations are the continuity equation and the London equation for 9. Here we outline 
the computation of £ . 

Since s C cS = n(—d t 6 — (V6) 2 /2m) + . . . one sees that the number density n is conjugate 
to — 9, and the energy density is given by 

£ = —nd t 9 + /in — £, (5.1) 

where dt9 is determined in terms of n, V8 and V^t by assuming a derivative expansion 
for dC/d{d t 6) = —n. From the effective Lagrangian at unitarity obtained in the previous 
Section 

C = c m 3 / 2 X 5 / 2 + Cl m 1 / 2 X- 1 /2(VX) 2 + c 2 m" 1 / 2 X 1 / 2 [(V 2 fl) 2 - ${d&0) 2 ] , (5.2) 
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one finds the energy density 

S = nV cxt + ~ n b/i + n 



3-2 2 / 3 (W) 2 8ci (Vn) 2 

cxt "I TFT- n + n ~ " ~ " — 

^/3 ( // ii m 2m 45com n 

~ "'- "^[(Vey-mW)*]- (5-3) 



Note that the next-to-leading order contribution to £ is exactly — £nlo if is replaced by 

o2/3 2/3 fVe) 2 

which is the solution of dC^o/ d(d t 9) = —n. 

The variation of the energy functional H[n, 9] = J d 3 xS with respect to 9 yields the 
continuity equation for the particle number 

SH „ dS _ / V9 \ 

66 (x) o(v6(x)) \ m J 

and the hydrodynamic equation for 9 is given by 

9 '* = -^)- < 5 - 6 > 

These equations describe the irrotational hydrodynamics of the superfluid at zero tem- 
perature. The equilibrium state in the presence of an external potential corresponds to 
9 = —[J,ot + cons, which gives a stationary density profile and zero superfluid velocity, 
V# = 0. The equilibrium particle density satisfies the condition 

22/3 2/3 3 2c i V 2 (V^) 

where //q is the chemical potential. Thus the leading behavior of n, which is obtained by 
dropping all the q, is given by a Thomas- Fermi approximation, whereas c\ determines the 
quantum kinetic energy correction and, as shown in AppendixEl the last term is proportional 
to the square of the traceless strain rate tensor. 

In writing Eq. (15. 2D we have used the results of the previous Section, namely = 
and /a//2 = —3. It is worth mentioning that, had we found a nonvanishing value for 
/4, the effective Lagrangian (15. 2p would have contained an additional term proportional to 
m 1//2 X 1//2 V 2 \4 xt , giving rise to a term proportional to 7?r 2 / 3 V 2 Kxt i n (15. 71) . Such a term 
would be incompatible with the boundary condition for n at infinity, due to the negative 
power of n. If we consider for instance the isotropic harmonic trap, where V 2 V ext is a 
constant, it is not possible to satisfy the asymptotics n — ► as r — ► oo. But, as mentioned 
above and shown explicitly in Appendix B, the condition f± = is an exact consequence of 
the conformal invariance of the field equations and the problem with the boundary conditions 
does not arise. 

Very recently, Rupak and Schafer j25| have derived an energy density functional using an 
epsilon expansion around d = 4 — e spatial dimensions. They find an expression depending 
only on n given by 

n 5/3 (Vn) 2 

£Rupak,Schaffer = ^xt™ + 1-364 + 0.032^ + O^Tl), (5.8) 

m run 
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which follows from their value of c\ — —0.0209 up to 0(e hie), together with £ ~ 0.4754 (or 
Co — 0.1166) in the same approximation j28|. Our result for £ is obtained by inserting the 
mean field results c = 2 5 / 2 /(15tt 2 £ 3 / 2 ) = 0.0842, Cl /c = -35/192 and c 2 = -0.0035 into 
Eq. (15.31) . We find that the energy density takes the form 

£ = y ext n + 1.6956 + n- ^ + 0.0324^ '— 

m 2m mn 

+0.0059— ((V 2 #) 2 - 3(d i d j 9) 2 ) . (5.9) 

It is remarkable the agreement between the terms corresponding to the quantum kinetic 
energy or quantum pressure, although they have been computed using two very different 
approaches. Although the values for the individual coefficients c and c\ show differences of 
the order of 30%, surprisingly these differences cancel in the ratio ci/c giving the quantum 
pressure. 

It is interesting to consider the expression for the energy density in the BEC limit which 
follows from the effective Lagrangian. It can be written in the form 

£bec = nV ext + ^n 2 + n^- + + + 2vra 3 n) V V ext , (5.10) 

2m 2m 32 mn 9bira 

where the most important derivative term is determined by the expression for fi(X) given 
in Eq. f!4.35bj) . If we make the replacements n — > 2ub, rn — > m#/2 and a —>■ Ob/2, this 



expression becomes 

S BEC = n M + + , B <™ + ^(V«' + . . . , (5.11) 

and one recovers the correct quantum pressure in the last term. Thus, the derivative part 
of the effective Lagrangian in the BEC limit fits in with the hydrodynamic description of a 
superfluid at zero temperature with bosonic constituents of mass 2m. It is worth pointing 
out the excellent numerical agreement between the coefficient 1/32 ~ 0.031 in the BEC limit 
and the corresponding coefficient 7/216 ~ 0.032 in the unitary limit. This is probably a 
reflection of the fact that, in the d — 4 limit about which the epsilon expansion is taken, the 
fermion-fermion scattering amplitude is saturated by the propagator of a boson with mass 
2m. 

Finally, in order to check the quality of the above NLO results in the BCS limit, a 
comparison with the predictions from the approach of Furnstahl et al. [29] would be very 
insightful. 



VI. CONCLUSION 

In this paper we have considered the derivative expansion of the effective action for the 
Goldstone field of a nonrelativistic superfluid Fermi gas at zero temperature in the whole 



BCS-BEC crossover. Based on the pioneering analysis of symmetries in Ref. |23j, we have 



shown that the NLO action can be given in terms of the four coefficient functions in 

c = p{x) + h(x)(yxy + ; 2 (x)(v 2 ^) 2 + / 3 pO(W) 2 + u(x)v 2 v ext , (6.1) 
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and we have given the precise relationships — see Eqs. (13.201) — between these functions 
and the derivatives of the response functions. It turns out that the computation of the 
NLO action relies also on the ratio of a pair of three-point functions, whose value we have 
determined by a calculation at the one-loop level. 

An important step towards Eqs. (13.201) has been the reduction of the initial La- 
grangian (13.31) to its final form ( 16. lft without higher-order time derivatives. In this regard, 
we have shown how the proper use of the LO field equations to eliminate higher-order time 
derivatives insures the consistency of the reduction process. As novel consequences of this 
critical analysis, we find the presence of a term proportional to ipidjO) 2 in (16. ip and the 
form of this action in the unitary limit 

C = cm V2 X 5 ' 2 + Cl m l ' 2 X- l / 2 {VXf + c 2 m- l ' 2 X 1 / 2 [(V 2 6) 2 - 3(d l d J 6) 2 ] , (6.2) 

which, as we have shown, is dictated by conformal invariance. In particular, conformal 
invariance of the field equation prevents the existence of an explicit coupling to the exter- 
nal field. 4 This Lagrangian determines uniquely the form of an energy density functional 
depending on the particle density n and the Goldstone mode 9, 

3-2 2 / 3 5/3 (V9) 2 8 Cl (Vn) 2 

£ = nV cxt H n/o n 5/3 + n- 



55/3 c 2/ 3 m 2m 45com n 



2 1 / 3 c 2 



1/3 



[{V 2 9f - mdjO) 2 ] • (6.3) 



It is worth mentioning that these aspects of our work do not rely on the specific approxima- 
tion used to compute the 2-point functions. Rather, they arise as an application of effective 
field theory ideas and techniques. 

By resorting to the one-loop approximation and taking into account the fi dependence 
of A in the mean-field approximation, we have also obtained analytic, closed expressions 
in terms of complete elliptic integrals for the coefficient functions of the NLO effective 
Lagrangian in the whole BCS-BEC crossover. Having closed expressions for these functions, 
rather than purely numerical results, makes the analytic study of different limits feasible. 
We have obtained series expansions near the unitary limit, and in the extreme BCS and 
BEC regimes. In particular, we have determined the mean-field values for the coefficients c\ 
and C2 at unitarity, and thus we have explicitly checked the form of the Lagrangian (16.21) . 
In this regard, it is interesting how the one-loop result found for the ratio of three-point 
functions combines with the result for g 4 to yield the conformally invariant combination 
X l / 2 [(V 2 d) 2 - ^(didjd) 2 )} —see Eqs. flSHD and (ET201) . Furthermore, the good agreement in 
the extreme regimes with known results from other approaches suggests that our mean field 
approximation can be taken as a reliable, first qualitative estimate for the coefficients of the 
effective theory. 
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APPENDIX A: ON THE REDUCTION OF THE INTEGRALS TO CANONICAL 
FORMS 

Here we outline the method used to avaluate all the momentum integrals in this paper 
in terms of elliptic integrals of the first and second kind. For the sake of concreteness, we 
consider the integral for Xm(0) 

Xnl(O) = , (Al) 

and assume \x > 0. The change of variables x = k 2 /2mfj, brings the integral into the following 
form 

m 3 / 2 A f°° , R(x) 
where y 2 = x[{x — l) 2 + a 2 } and R(x) is the rational function 

«M - • < A3 > 

Now, the integral of the quotient of a rational function by the square root of a cubic 
or a quartic polynomial is, by definition, an elliptic integral. But we still have to reduce 
this integral to a combination of canonical forms. The fact that the integrand has no other 
singularities than the branch points of y implies that only elliptic integrals of the first (F) 
and second kind (E) can be involved. These are given by [30] 



F(<p\n) 



sin ip 



dt 



o V(l-t2)(i_nt2) 



E,sVn = i <ll^=- (A4) 



The change of variables 



t = x-^ + a 2 fA5 , 
x + VTTa 2 ' 



brings the integration interval to (—1,1) and y into the canonical form 

dx \[2 dt 



V V-l + Vl + a 2 ^(l-t 2 )(l + 1 t 2 

where 



(A6) 



7 = 1 a • A7 

or 



Finally, after using standard "reduction formulae" 3(| we obtain 



™3/2 A i 

Xnl(0) = - , 2 , 1/2 ° ; ==if(-7) , (A8) 



a 2 
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where K(—j) = F(jr/2\ — 7) is the complete elliptic integral of the first kind. 

The evaluation of the other integrals proceeds along the same lines. In all cases we have 
an odd power of Ef. = fi[(x — l) 2 + a 2 } that combines with yfx from the measure d 3 k to 
give an integrand of the form R(x)/y, where only the rational function R(x) changes from 
case to case. As the only singularities are the branch points of y — there are no additional 
poles — only complete elliptic integrals of the first and second kind can result. 



APPENDIX B: CONFORMAL INVARIANCE OF THE NLO LAGRANGIAN 

In this Appendix we investigate the constraints imposed by conformal invariance on the 
form of the NLO Lagrangian in the unitary limit. As shown in Refs. (23, 24j, in that limit the 
original fermionic action 5 is invariant under a special set of time- dependent transformations. 
Actually, Son and Wingate [23[ have discussed the invariance properties of the action when 
the system is put in a curved manifold and an external gauge field. In order to adapt their 
transformations to our case, where the metric is euclidean and the gauge field is set to 
zero, we have to use a "gauge fixed" version wich includes, besides the purely conformal 
transformation, "compensating" coordinate and gauge transformations. Its infinitesimal 
form is given by 

Si> = /3"(0^^ - ^x • W - WW ~ (Bl) 

SV cxt = -/T(t)^ - ^-x • W e xt - P(t)dtV ext - P\t)V ext , (B2) 

where (3{t) is an arbitrary function of t. In the notation of Ref. (23|, this transformation can 
be seen as a combination of general, gauge and conformal transformations with parameters 
£ fc = P'(t)x k /2, a = m(3"{t)x 2 /4 and (3(t), respectively. With this choice, the euclidean 
metric and the zero external gauge field are untouched. We assume that the variation of 
the chemical potential is given by 5fi = —j3'(t)fi, which assigns to V ext all the change in the 
variation 5(V ext — yu) under a gauge transformation. In the effective theory the relevant field 
is the Goldstone mode 6, which transforms inhomogeneously according to 6 

se = jT(t)^--^-x-ve-p(t)dte. (B3) 

A "scale transformation" is a particularly simple conformal transformation where /3(t) is 
a linear function. As shown in Ref. [23[, scale invariance alone determines the form of the 
functions fi(X) in the NLO Lagrangian (14.251) 

£nlo = cm^X-^iVX) 2 + m-^X 1 ' 2 [c 2 {V 2 9) 2 + c 3 (d t d 3 6) 2 + c 4 mVV ext ] (B4) 

The change under a general conformal transformation is then given by 

5£ NLO = (3c 2 + c 3 )m 1 / 2 X 1 / 2 V 2 9(3"(t) - ^cm^X 1 / 2 (3"\t) 

-d t (£nlo/9(*)) - • (x£ NLO p'(t)) . (B5) 



5 By this we mean the fermion action before the Hubbard-Stratanovich transformation leading to Eq. (|2.ip 
is applied. 



The Goldstone field here corresponds to fit — #son,wingatc of Ref. 23 1 
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Note that, as expected, the action is automatically invariant under scale transformations, 
for which (3"{t) = (3"'{t) = 0. Invariance under general conformal transformations requires 

c 3 = -3c 2 , c 4 = , (B6) 

while the value of C\ is unrestricted. This constrains the NLO Lagrangian to the form given 
by Eq. flQSP . 

Note that we could integrate the variation of the c 4 term in (IB5I) by parts to get a 
contribution proportional to f3"(t). Using the LO field equations to eliminate the time 
derivative d t X would then give the following result 

5£nlo = (3c 2 + c 3 + c 4 )m 1 / 2 X 1 / 2 V 2 6 f3"(t) (B7) 

up to total derivatives. This suggests that the action is conformally invariant as long as 

c 3 = -3c 2 - c 4 , (B8) 

which is a weaker constraint than ( 1B6I) . But here, as was the case with the elimination of 
the 'g^{X)X t V 2 9 term at the end of Section [TTTl the use of the LO field equation is not fully 
legitimate. The reason is that the elimination of the time derivative dtX in favour of V 2 6* 
would involve dropping a term proportional to the LO field equation, namely 

^rn*X-^nt) (X t + gv 2 *) . 

It is easy to check that the Euler-Lagrange equation arising from this term gives a contri- 
bution to the linearized field equation proportional to 

3md t [/?'"(V ext + d t 6)\ + 10^"V 2 6 - 4/i/TvVext, (B9) 

which can not be neglected. In other words, the field equations derived from a Lagrangian 
subject only to the weaker constraints flB8j) will not be invariant under conformal transfor- 
mations. Thus, the coefficient of (3'"{t) in Eq. ( 1B5I) must be zero irrespective of the values 
of c 2 and c 3 . This fact, which seems to have been overlooked by the authors of Ref. [23~l ], 
would explain the discrepancies with their results. 7 

We end this Appendix by noting that the the conformally invariant combination 
(V 2 6*) 2 — 3(didj8) 2 can be written as the square of an I = 2 irreducible rank-two tensor 
of conformal dimension one 

(V 2 #) 2 - ^didjOf = -3 (didfi - i<%V 2 ^) . (BIO) 

This allows a natural interpretation in terms of superfluid hydrodynamics. Writing the RHS 
of (IBlOp in terms of the superfluid velocity v s = m _1 V6 l shows that the new invariant is 
proportional to the square of the traceless part of the strain rate tensor, also known as the 
shear rate tensor 

-(diV s j + djV 8)i ) - - 5ijV ■ v s . (Bll) 



7 See also comments at the end of Section IIII1 
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APPENDIX C: DERIVATIVES OF THE ONE-LOOP RESPONSE FUNCTIONS 



Here we give the derivatives of the response functions used in the computation of the 
NLO action at the one-loop level. The results are 

OXnn f d 3 k Al 



dto 2 " J (2tt) 3 4£| 

m 3/2 y I 

= 24^^ a 2 + <* V - SignM + - E ^ 

+ (3 + 3a 2 - sign^Vl + a 2 ) #(-7)] , (CI) 

dxnn f d 3 k / 3Ag& A 2 (Al-4e k )k 2 \ 
dq 2 ' J (2vr) 3 \8mEl 2Am 2 E 7 k ) 

m i/2 x 



247r 2 |/i| 1 /2 i + a 2V si § n ^) + Vl + a 2 [sign(/i)£(- 7 ) 
+v / TT^ 2 ^(-7)] , (C2) 



dxm f d 3 k A £ k 



/• d A k 

J W) 



du 2 ' J (2tt) 3 AEl 
m 3 / 2 A 1 



-sign(/i) + Vl + a 2 [— sign(//)£'(— 7) 



247T 2 |/i| 5 /2 a 2 + a 4 

-v / TT^ 2 ^(-7)] , (C3) 

Sxni /• gfe / A (2el-A 2 ) A q(3A 2 -24 2 )A; 2 \ 
<9g 2 7 (2vr) 3 V 8mE| 24m 2 £ fe 7 / 

= 24^|^ a 2 I a« \f~^ + v/f + ^ 

- + a 2 + sign(/i)v / TT^ 2 ) ^(-7)] , (C4) 
9xn f d 3 k H 



du 2 J (27r) 3 4E| 

m 3/2 l 



~ 2A, 2 W 2 a 2 + * V"^) + «" 2 " 

+ ^3 + 3a 2 + sign(/i)v / rT^ 2 )^(-7)] , (C5) 

9xn = /■ rf 3 fe / &(ff-2Ag) 5Agffi» \ 
dq 2 J (2vr) 3 V 8mE| 24m 2 E 7 k J 



m 



1/2 



: ^/-sign(//) + Vl + a 2 [(4 + a 2 )sign(/i)£(-7) 



727r 2 |/i| 1 /2 a 2 + a 4 

+ (6 sign(/i)(l + a 2 ) + (10 + 7a 2 )Vl + a 2 ) ^(-7)] , (C6) 
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d X L _ f d 3 k A 2 k 4 



dq 2 ' J (2tt) 3 20m 4 £f 

1 1/2 - 



+ (3 + 3a 2 - sign(/i)v / TT^) K(-y)] , (C7) 
0, (C8) 



9xl 
duj 2 

d Xjn f d*k A 2 k 2 



duj J (2tt) 3 Ylm 2 E\ 

m l/2 J 



sign(/x) + vl + a 2 [-4 sign {\i)E (-7) 



367r 2 1 I x / 2 a 2 

+ (3 sign(/i) - v/TT^ 2 ) AT(— -y)] , (C9) 



dxji f d 3 k A £ k k 



2 



duj J (2tt) 3 Ylm 2 E\ 

m 1 / 2 A 1 



12^1^3/2 a 2 V Sign(/i) + V/IT ^ hjB( - 7) + ^ ( - 7)] • ' ( :1 ° ^ 

The derivatives of the transverse current response are in the same approximation 

9xt 1 



dq 2 3 <9g 2 



(Cll) 



du 2 



0. (C12) 



APPENDIX D: ONE-LOOP THREE-POINT FUNCTIONS AND THE RATIO 

53/52- 

In this Appendix we show how the computation of a pair of three-point functions in 
the one-loop approximation leads to the result g 3 (X)/g2(X) = 2. In order to identify the 
coefficients of (V 2 #) 2 and (didj6) 2 , we first write down all the third-order terms proportional 
to Y66 in the effective Lagrangian (13.31) 

= —g^Ve ■ VYV 2 9 - g'MY(V 2 ef - g'MYid^ef, (Dl) 

where Y = \x — X. This cubic Lagrangian gives rise to a coupling Y99 with vertex propor- 
tional to 

-gM [qi ■ Q2{ql + ql) + H4] - ^MqUI ~ ^M{0.\ ' <? 2 ) 2 , (D2) 
m 

where q\ and qi are the momenta of the pair of Goldstone fields. Note that the coefficient 
of the first term is known in terms of two-point functions 

95 = „ 6 -6 3 = -^ + ^_^i. (D3) 
ouo 1 + gxw ouj 
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Next, we compute the Fourier transforms of the one- loop three-point functions that induce 
the effective vertex f]D2|) . These are given by (T |\l/" l "(x)r 1 ' 3 \l/(x) V • j p (y)V • j p (z)}}, where 
j p is the paramagnetic fermionic current that couples to V# in Eq. (I2.3j) . The two triangle 
diagrams that contribute read 

rfliQuQi) = \ tr [g(K)g(K + Q 1 )g(K + Q 1 + Q 2 )r 1 > 3 ] 
Jk 



x [k-q 1 + ^\(k-q 2 + q 1 -q 2 + f\^ + (Q 1 ^Q 2 ). (D4) 



For our purposes, it is sufficient to obtain the 0(q 4 ) contribution to T^ 3 ] at zero frequencies. 
A lengthly computation yields 



d 3 k A 2 k 2 {k 2 £ k -2mEl] 



,2 „2 



+ J (2tt)3 Um^El qiq2 

f d 3 k A 2 k 2 (k%-2mE 2 ) Q? + 0(a5) (m) 

+ J J^y e^Ei {qi q2) +0(g) ' (D5) 



r (3) = 



^T2^Et^- q2 ^ +q2) + 2q ^ 



+ J (2tt)3 Q0m*El 9l?2 

+ J J^y so^ej {qi ■ q2) +0{q) - (D6) 

We observe that the coefficients of (q 2 + q 2 ) + 2qfq 2 in and are given respectively 
by —m~ 1 dxjn/du} and —m^dxji/duj (see Eqs. (IC9j) and fIClOl) ). This agrees with Eq. (1D3I) . 
after the replacement a — > gxnY / (2+ gXu) — from the gap equation — is made in the term of 
the effective action proportional to vf^aOO. Simple inspection of the coefficients of (qi ■ q 2 ) 2 
and q\q\ shows that g' 3 (f^)/g 2 (^) = 2, and therefore g 3 (X)/g 2 (X) = 2, after dropping an 
irrelevant integration constant. 

The result g 3 {X) / g 2 (X) = 2 is strikingly simple, and one might wonder whether this is 
just a peculiarity of the one-loop aproximation used here. But, unlike the ratio f 3 J f 2 = —3 
found in Section [TTT1 the value g 3 (X) / g 2 (X) = 2 does not seem to be the consequence of any 
obvious symmetry. 
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